On The Geometrical Description of Dynamical Stability 
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A general parametrization requirement for geometrization of dynamical stability is presented. 
We show that the non-physical behaviors appear in lower dimensional systems when non-affine 
parametrization of arc length with time is used. We compare the two widely used Jacobi and 
Eisenhart metrics as archetypes for (non)a£Bne parametrization. We numerically investigate this in 
the context of the two-centered Morse potential. The relevance of parametric resonance as a source 
of instability in two dimensional systems is resolved. 



PACS numbers: 



The general time-independent Hamiltonian 
H{q,p) = ^a'^q)p^PJ+V{q) 



(1) 



is the basic representation of many problems and appli- 
cations in the physical sciences ranging from molecular 
physics to astrophysics. Naturally different approaches 
for fundamental understanding of such systems and their 
dynamical (in) stabilities have been the subject of vigor- 
ous research P, Due to its non-perturbative nature, 
geometrization of dynamics has been proposed as a novel 
approach , which gives more physically meaningful 
results for the stability when the traditional Lyapunov 
exponent does not.[5| 

The motion of a mechanical system is entirely deter- 
mined by the principle of least action Q. Geometriza- 
tion of dynamics represents the motion as a geodesic on 
a manifold with a suitable metric g — 
For a congruence of geodesies 7'(t, s), the vector field 
of geodesic spread is defined as 



dYiT,s) 



dr 



(2) 



where s is the arc length and r parameterizes the dif- 
ferent geodesies. £_{s) is a measure of the separation of 
nearby geodesies and its evolution in time or arc-length 
has been considered as an indicator of (in)stability. This 
approach can be utilized to tackle a broad range of ap- 
plications in non-integrable systems and can give new 
insights to the on-set of chaos. evolves according to 
Jacobi-Levi-Civita (JLC) equation^], which in local co- 
ordinates reads as 



~d^ 



dq^kdq" 



ds ds 



= 



(3) 



where ^ and i?*j.^stand respectively for the covariant 
derivative along geodesies and the Riemann-Christoffel 
curvature tensor. Among other possible choices, two met- 
rics that have been extensively used are: the Jacobi or 



the kinetic energy metric (gj) on the admissible configu- 
ration space Mj = {x e M : V{x) < E} 



{gj),j=2[E-V{q)]a,,{q) 



(4) 



and the Eisenhart metric 8]. The latter is defined on the 
enlarged (M x 3?^, gs) configuration space-time x" — 
{q° = t, g*, with the arc-length 



ds' 



2V{q){dq'y + aijdq'dq^ + 2dq''dq 



(5) 



In the Eisenhart 's metric physical motions satisfy an 
affine parametrization ds^ = K^dt^ {k is a real constant). 
The time evolution of the extra coordinate can formally 
be integrated 



-,N+1 



—t + Co~ / Cdt' 
2 Jo 



(6) 



Equations of motion for ^ using Eisenhart metric are 
equivalent to those of the tangent dynamics, thus show- 
ing its geometrical origin Solutions of JLC equa- 
tions are quite complicated and computationally expen- 
sive to obtain; they reqiiire the simultaneous computa- 
tion of (7*(s) and ^*(s) [10|. For this reason numerical 
integration of the exact equations ^ has been mostly 
restricted to low dimensional systems (i.e. two degrees 
of freedom)!^, 11 1 and the equivalence of results using 
both metrics has been argued by simulations of particu- 
lar systems]^, [l3|. In the case of many degrees of free- 
dom, an approximate form of JLC equation for ||^|| has 
been obtained 10]; this scalar function supposedly gives 
sufficient information about stability. Moreover, in the 
thermodynamic limit, under statistical considerations a 
version of this treatment has been proposed whose main 
goal is to quantify the stability based only on static ge- 
ometrical properties [H, [13. This by -passes the need 
of computing trajectories for a representative set of ini- 
tial conditions and gives rise to a potentially powerful 
tool for finding order-to-chaos transitions. The main ar- 
gument in favor of using Eisenhart metric up to now is 
its computational simplicity. On the other hand in the 
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framework of Jacobi metric, the scalar, Ricci's and sec- 
tional curvatures have been shown to be positive in many 
systems of interest [13, This led to the conclusion 
that negative curvature is not the fundamental source of 
instability and to the conjecture that the oscillating cur- 
vature leads the systems to chaos through parametric res- 
onance. Moreover, Jacobi metric has been considered to 
be a preferred choice due to its mathematical properties, 
namely its compactness and Riemannian character [l^ . 
It has been used for stability analysis of a system of self- 
gravitating particles [lOf. Interestingly it was observed 
that, for the same energy per degree of freedom, increas- 
ing the number of particles leads to less divergence in the 
geodesic spread namely suppresses chaos. These result 
of simulations were assumed to be due to the mathemati- 
cal approximations [lol . 11 1 that were made for obtaining 
the dynamics of ||^||. Instead, one expects the probabil- 
ity of chaotic dynamics to increase with the number of 
non-separable degrees of freedom N. 

In the present letter we wish to show the fundamentally 
different results that can emerge from the application of 
the different metrics mentioned above. We show that in 
the low dimensional systems non-affine parametrization 
naturally leads to unstable solutions that should neither 
be confused with parametric resonance nor be interpreted 
as physical instabilities. We will discuss conditions under 
which stability analysis can be performed using either of 
the two metrics. 

A manifold has a well-defined tangent space every- 
where. Boundaries of accessible configuration space dMj 
in Jacobi metric correspond to points where the kinetic 
energy, T, is zero (i.e. E = V{q)). The physical picture 
is that the geodesies reflect from dMj. Moreover, the 
reflections have to be perpendicular to dMj otherwise 
there would be a net momentum parallel to dMj, which 
would contradict the definition of dMj (T = 0) [ll, [3l . 
The curvature tensor in terms of Christoffel symbols is: 



R 



ilk 



jk 



dr 



r 



jk 



jk n 



,k + 9kr, 



9jk,m} 



(7) 

(8) 



Here, gij is given by Eq.Q; thus Fjj, contain powers of 
[E — ^(9)]^^ that become infinite on dMj. Clearly this 
corresponds to cusps in the curvature, where the tan- 
gent space is not defined. Moreover there are station- 
ary geodesies at the boundaries, which in general are not 
equilibrium points. Therefore, these trajectories are not 
physically meaningful on dMj, in addition to the usage 
of manifolds not being allowed. However, in the Eisen- 
hart's metric all curvatures are finite (Ricci, sectional, 
and scalar). From a physical point of view the stability 
is determined by comparing the difference of trajectories 
in time, say from the value of Lyapunov exponent. 

It is evident from (Eql2|) that ^ by definition com- 
pares the geodesies at the same arc length. When 



parametrization is affine, as in Eisenhart's metric, com- 
paring geodesies at equal arc length corresponds to com- 
parison of trajectories at the same time. Therefore, 
the stability analysis is in a one-to-one correspondence 
between the divergence of trajectories in time and the 
geodesic spread. 

A natural question is where and under what condi- 
tions does the Jacobi metric give similar results for sta- 
bility analysis to those given by the Eisenhart metric? 
A first criterion would be to avoid singularities in the 
curvatures. A conserved nonzero angular momentum for 
bounded systems warranties that all momentums are not 
simultaneously zero, which puts a lower bound on the ki- 
netic energy. Another favorable condition could be a high 
number of degrees of freedom, as this generally reduces 
the oscillations of the kinetic energy. By the arguments 
given above, the usage of Jacobi metric method for sys- 
tems with few degrees of freedom and zero total angular 
momentum should therefore be avoided. A detailed anal- 
ysis with respect to the dimensionality and coherence of 
the system will be addressed in a sequel. Here we focus 
on a two-dimensional problem. 




FIG. 1; Poincare maps for AE — 7.5 

In order to make our analysis more illustrative we ap- 
ply this formalism to a two-center Morse potential with 
center-to-center distance r^ = 3.0271A, and Morse pa- 
rameters D = 40.75 cm"\ a = 1.56A~\ and J= 4.361 
0. The corresponding Hamiltonian with zero total an- 
gular momentum is 



H 



p2 



1 



1 



2^ii?2 2^i2rl 



Pi + V{R,9) (9) 



where, R and 9 are polar coordinates; the polar axis runs 
across the two centers, and the reduced masses are taken 
as fii = 18 amu, /X2 = 64 amu[9l|. In the figures, we 
have expressed all relevant quantities in units that are 
powers of cm (i.e. energy, linear momentum and A are 
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in cm ^, t and distances are in cm, thus angular mo- 
mentum becomes dimensionless) . In Fig|l] we show a 
Poincare Surface Section (PSS) for AE ^ E - 2D ^ 7.5 
cm^^ where it is predominated by stable limit-cycles, 
substituting for s — > s(t) We plot the Lyapunov expo- 
nent A = hmj^oo J log( 111(0)11 ) versus the energy using 
the Eisenhart metric (Figl2|). A is a function of energy 



15 



A Energy cm ^ 



FIG. 2: Lyapunov values versus Energy using Eisenhart's 
metric 



and the initial conditions; we choose a set of initial con- 
ditions in the semi-stable subspace for AE = 7.5 cm~^ 
in agreement with Fig[TJ As one expects the lower en- 
ergy values correspond to stable solutions whereas the 
higher energies, AE > 4.3 cm~^, produce positive val- 
ues of A and therefore, the solutions exponentially di- 
verge. In order to ensure a proper convergence to the 
Lyapunov exponent we show in Figl3] the evolution of 
Y(t) — J log( [1110)11 ) for different energies ranging from 
very stable {AE = 0.5 cm^^) to unstable (AE = 15.5 
cm~^) regimes. It is important to note that the stable 
dynamics demands A < 0. The same is shown in Fig H] for 
a very stable regime {AE = 2.5 cm~^) using Jacobi met- 
ric. As can be seen, Y converges to about 130 cm~^, in 
contradiction to the PSS and to the results obtained from 
Eisenhart metric (see Figl2])! We also plotted in FigElthe 
behavior of the Y function in a small time window for 
the same initial conditions and energy as in FigU) It is 
clear that the nonphysical limit for the convergence of the 
Y-function in FigH] is a consequence of the "kicks" pro- 
duced by the oscillations in the kinetic energy. Although 
we have chosen as an example a particular system, the 
appearance of singularities in the curvature and its oscil- 
lations are general consequences of the Jacobi metric. In 
this case, the curvature is actually positive and therefore 
the exponential growth of ^ is due to the dominant quasi- 
periodic "kicks" produced by the kinetic energy values 
close to zero as is evident in (equations [31 and [7]) . 
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FIG. 3: Lyapunov convergence for different energies using 
Eisenhart metric 
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FIG. 4: Kinetic Energy and Lyapunov convergence for 
AE—2.5 cm~^ using Jacobi metric 



Before concluding we make some final remarks to put 
our work into a larger context of recent developments. 
Though still mathematically unproven we have shown a 
counter-example against the reliability of the Jacobi met- 
ric for geometrization of dynamical stability. The reason 
this has gone unnoticed in the past is most probably due 
to the fact that the formalisms were applied in regimes, 
where i) Jacobi metric could be approximated as afRne, 
i.e. small changes of the kinetic energy, as can be the 
case for large number of degrees of freedom 12,13,17, 
ii) the particular conditions, under which comparison be- 
tween tangent dynamics and the Jacobi metric were done, 
were already in the chaotic regime or at the onset of 
chaos; thus the discrepancy between Eisenhart and Ja- 
cobi metrics could not be clearly established 13, 11 1 iii) 
calculations were done mostly using the Eisenhart met- 
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FIG. 5; Zooming in on Fig|4l Kinetic Energy and the corre- 
sponding "kicks" in ^ for A_E=2.5 cm~^ 

ric and assumed to hold for Jacobi as well (please see 
references in Furthermore, parametric resonance 

is a very particular mechanism for driving a system to 
instability (please see Landau and Lifshitz [6]), which 
in the past has not been firmly investigated as a source 
of chaos within the Jacobi framework. Here we show 
how semi-periodic kicks produced by small values of T, 
which is a side-effect of the usage of the Jacobi metric as 
opposed to intrinsic dynamical property of the system, 
may be misinterpreted as actual physical parametric res- 
onance. To our knowledge Eisenhart's metric remains to 
be a very successful method for geometrization of general 
holonomic conservative systems in classical dynamics. In 
addition geometrizatoin is much easier to work with using 
Eisenhart's metric. 

In this letter we have discussed an important require- 
ment for geometrization of stability analysis. We have 
explicitly shown that the time-evolution of the geodesic 
spread, in lower dimensional systems, is not a measure 
of stability when the parametrization of the arc length 
with time is non-afhne. Due to their importance, we 
put this discussion in the context of Jacobi and Eisen- 
hart metrics, which can be looked at as the archetypes 
of non- and afRne parametrizations respectively. The 
apparent parametric resonance phenomena seen in the 



frame-work of Jacobi geometrization discussed in earlier 
works is most probably due to the oscillations of the ki- 
netic energy regardless of the stability of the system. Fu- 
ture work entails a thorough investigation of stability for 
higher degrees of freedom systems as well as the explicit 
dependence of the stability on the degree of non-afhnity- 
the fluctuations in the kinetic energy. 
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